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We present preliminary results on the properties of relativistically broadened Fe K Q lines in a collection of more then 100 
Active Galactic Nuclei (AGN) observed by the XMM-Newton EPIC-pn camera. Our main conclusions can be summarized 
as follows: a) we detect broad lines in about 25% of the sample objects. This fraction increases to 42 ± 13% if we consider 
only objects with more than 10 4 counts in the hard (2-10 keV) band, and to 50 ± 32% for the small sub-sample (6 
objects) of type 1 Piccinotti AGN with optimal XMM-Newton exposure (at least 2 x 10 counts in the hard band); 
b) we find no significant difference in the detection rate of broad lines between obscured and unobscured AGN; c) the 
strongest relativistic profiles are measured in low-luminosity (Lx < W 43 erg s _1 ) AGN; d) Equivalent Widths (EWs) 
associated with relativistic profiles in stacked spectra are ^,150 eV for all luminosity classes; e) models of relativistically 
broadened iron line profiles (kyrline, Dovciak et al. 2005), which include full relativistic treatment of the accretion 
disk emission around a Kerr black hole in the strong gravity regime, yield an average disk inclination angle ~30°, and a 
radial dependence of the disk emissivity profile ~-3. The distribution of EW is very broad, with (log(EW)) = 2.4 and 
<7i og (EW) — 1-4. We estimate that an investment of about 1 Ms of XMM-Newton time would be required to put these 
results on a sound statistical basis. 

© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 

1 Introduction proach that has inspired the study presented in this paper. 

We aim at answering primarily the following questions: 

XMM-Newton is the most suitable mission to study the prop-1 

— how often are effects due to X-ray illumination of a rel- 
ativistic accretion disk detected in AGN? 

— what do the properties of relativistically broadened pro- 
files tell us about the accretion flow in the innermost 
regions around super-massive black holes? 

1.1 The sample 

The sample comprises 102 AGN observed as pointed nom- 
inal targets by XMM-Newton, and whose data were avail- 
able in the public archive as of June 2006. Out of them: 32 
are classified as Seyfert Is, 18 as Narrow Line Seyfert Is 
(NLSyls), 39 are Radio Quiet Quasars (RQQs), and 13 are 
Seyfert 2s. No attempt has been made to create a complete 
and unbiased sample. We have restricted our study to Compton 
thin obscured Seyfert 2s, whose nuclear emission is covered 
by a column density Nh < 10 22 5 cm~ 2 , because the de- 
termination of the K a iron line profile in spectra affected 
by stronger obscuration is uncertain. Most of the objects are 
local, with ~60% of them having a redshift < 0.1, and 98% 
of them < 0.5. 

2 Data reduction and analysis 



For each of the observations, data have been reduced from 
the Observation Data Files using SASv6.5 (Gabriel et al. 



erties of broad iron K Q emission lines, thanks to the un- 
precedented collecting area of its EPIC cameras. As this 
workshop has demonstrated, the observational properties of 
the Fe emission line complex have been one of the primary 
subjects of AGN studies with XMM-Newton (see Fabian, 
this volume, and references therein). Although most of the 
papers in this field have focused on the detailed properties 
of the iron line profile in individual bright AGN, there have 
appeared in the literature some attempts to estimate the sta- 
tistical importance of relativistic effects. Stacking spectra of 
a large sample of z ~ 1 AGN in the Lockman Hole, for in- 
stance, Streblyanska et al. (2005) unveiled a relativistically 
broadened and skewed profile in both type 1 and type 2 
objects. The Equivalent Width, EW, of such a feature is 
large (in the range 420-560 eV for unobscured, 280-460 eV 
for obscured AGN, depending on the profile models used). 
However, broadened iron K a lines with such a large EW are 
uncommon in the local universe (cf., e.g., the distribution of 
the broad Fe line EWs in a sample of 40 PG quasars studied 
by Jimenez-Bailon et al. 2005 in Fig.Q. 

After almost seven years of successful XMM-Newton 
operations, it's time to give a systematic look backward at 
the wealth of data contained in the archive. This is the ap- 
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Fig. 1 Distribution of the Fe K a line broad component 
EW in the sample of 40 PG quasars studied by Jimenez- 
Bailon et al. (2005). Shaded cells represent 90% confi- 
dence level upper limits. The measured EW in Lockman 
Hole unobscured AGN [after Tab. 1 in Streblyanska et al. 
(2005)] is marked by the hatched rectangle [420±|§ eV 
when the broad line is modeled with a Gaussian profile, as 
in Jimenez-Baildn et al. (2005)]. 



2003), and the most advanced calibration files available as 
of March 2006. We present in this paper data from the pn 
camera (Striider et al. 2001) only, given its larger effective 
area. Spectra were extracted using standard pattern event 
selection. The radius of the circular spectral extraction re- 
gions and the count rate background rejection thresholds 
were chosen in order to maximize the signal-to-noise in 



Fig. 2 EW of the relativistically broadened Fe K a line 
against the net hard {i.e. 2-10 keV) counts. Each data 
point corresponds to one XMM-Newton observation. Mul- 
tiple observations of the same objects are shown as dis- 
tinct points. Filled squares represent line detections, empty 
squares represent upper limits. Source excluded from the 
sample are labeled (see text for details). The three objects 
with EW upper limits <20 eV are: IC4329A (see, e.g. Gon- 
doin et al. 2001, Steenbrugge et al. 2005), Mkn 509 (Pounds 
et al. 2001), NGC7469 (Blustin et al. 2003). 



centroid energy 6.4 keV (Fel), 6.7 keV (Fexxv), and 6.96 keV| 

(FeXXVl), and B is a relativistically broadened profile. We 

have used for B model profiles as predicted around Schwarzschild| 

(a = 0; Fabian et al. 1989) and Kerr (a ~ 1; Laor 1991) 



each individual spectrum. Background spectra were extracted! black holes, as well a set of models (kyr 1 ine; Dovciak et 
from nearby circular areas in the same chip as the source, al . 2 005), which includes full relativistic treatment of the ac- 
except for observations in Small Window mode, where blank| cretion disk emission around a Kerr black hole in the strong 



field event list (Read & Ponman 2003) were used, after co- 
ordinate recasting and normalization to the quiescent back- 
ground level. For the time being, we have not merged differ- 
ent observations of the same target. Whenever multiple ob- 
servations are available, we have used that with the longest 
"good" exposure time (except in Fig. [2] see Sect. 3). 

All the spectra have been rebinned, in order to ensure 



gravity regime. The ky model family depends on the disk 
inclination angle tp, on the radial dependence of the emis- 
sivity (power-law index (3) and on the dimensionless black 
hole spin a. In order not to over-fit the data in moderate sta- 
tistical quality spectra, we have employed the (3 parameter 
only as a measurement of the strength of relativistic effects, 
while fixing the innermost and outermost radii of the line- 



that each background-subtracted spectral channel has atleast| emitting region to the innermost stable circular orbit for a 



25 counts, and that the instrumental energy resolution is not 
oversampled by a factor larger than 3. Spectra have been 
fit in the 2.5-15 keV energy range. We have employed the 
following baseline fitting model: 



x A[E~ T + C(r, </>, E) + Y? l=l G l + B((f> 1 0, a)](l) 



given a value, and to 400 gravitational radii, respectively. 

3 Detection of relativistic lines in individual 
spectra 



where dp/j is the photoelectric cross-section (Anders & Grevessi( Fig. |2l we show the EW of the Fe K Q relativistically 
1989), Nji is the column density (constrained not to be broadened component as a function of the total net counts 
lower than the contribution due to intervening gas in our in the 2-10 keV band. As expected, the overwhelming ma- 
Galaxy), C is a Compton-reflection component by a plane- jority of the measurements are aligned along a straight line 
parallel neutral slab (model pexrav in XSPEC; Magdziarz (in log-log space), which represents the sensitivity limit of 
& Zdziarski 1995), Gi are Gaussian unresolved profiles with| our sample. Three poorly exposed (< 10 000 hard counts) 
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Fig. 3 Distribution of the net counts in the hard X-ray 
band for the AGN of our sample. Filled cells correspond to 
detections of a broad iron K Q profile; shaded cells to non- 
detections. 



AGN type 


Well exposed 


Underexposed 




(%) 


(%) 


Seyfert Is 


50 ±15 


< 21 


RQQs 


60 ±48 


17 ±8 


NLSyl s 


< 30 


20 ± 14 


Seyfert 2s 


33 ±24 


< 30 



Table 1 Fraction of broad iron K a detections in our AGN 
sample. Errors are at the la level 
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Fig. 4 Distribution of upper limits on the EW of the broad 
iron K a line in our sample. 



Parameter 





a 


N obj 


P 


-2.6 


1.8 


22 




33 


21 


21 


log(£W) 


2.4 


1.4 


20 



Table 2 Mean and standard deviations on j3, (j>, and EW 
in our sample. N bj is the number of objects on which a 
statistically meaningful constraints on the parameter can be 
set, and which have been used to construct the normalized 
distributions in Fig. [5] 



objects are well above this correlation, and exhibit EW by 
an order of magnitude larger than the average of the sam- 
ple at their net count level. They are: IRAS 13224-3809, 
PG141 1+442, and 1H 0707-495. In the last the presence of 
a EW > 1 keV relativistically broadened Fe K a line has 
been challenged (Gallo et al. 2004, Fabian et al. 2004). We 
therefore consider these detections as "suspicious", and re- 
moved these objects from the sample. For similar reason, 
we exclude from the sample PG121 1+143 (Pounds et al. 
2003). None of the sample results presented in this paper is 
significantly affected by this choice. 

The distribution of the source counts is shown in Fig. [3] 
In these histograms, filled cells represent objects where a 
broad Fe K a iron line is detected, shaded cells represent 
non-detections (the distribution of the EW upper limits is 
shown in Fig. |4}. The fraction of detections is 42 ± 13% 
and 16 ± 6% for "well-exposed", and "under-exposed" ob- 
jects, respectively, where we have arbitrarily set a threshold 
of 10000 net counts to distinguish the two classes. The de- 
tection fraction is not significantly dependent on the AGN 
type (Tab.[0; however, the small number statistics prevents 
strong statement from being drawn. 

In Fig. |5] normalized distribution functions are shown 
for the best-fit parameters on which the kyrline model 



depends. These distributions are generated combining one 
Gaussian distribution for each individual measurement, whosel 
mean (standard deviation) is the best-fit value (statistical er- 
ror ). In Tab.EJthe mean (()) and standard deviation (cr) of the 
combined distributions are listed, together with the number 
of measurements used to construct the normalized distribu- 
tion. In 5 objects a constraints on the black hole spin can be 
derived ( (a) = 0.6, a a = 0.3). In two of them (HE1029- 
1401, a = 0.85 ± 0.04; NGC 3227, a < 0.18) the measure- 
ment is formally inconsistent with a maximally spinning 
black hole; in the other three (MCG-6-30-15, NGC 3516, 
NGC 4395) only a lower limit on a can be derived. 

Several authors have discussed the effect that the spec- 
tral curvature produced by high-density ionized absorbers 
may have on the detectability of broad Fe K Q profiles (Reeves| 
et al. 2004, Pounds et al. 2003, Turner et al. 2005; see as 
well Diaz-Trigo's and Reeves' papers in this volume). If not 
properly accounted for by "warm absorbers" models, this 
curvature may mimic the red wings of a relativistic broad- 
ened profile in the 3-6 keV energy range. Approximating 
an ionized absorber with a cold absorber - as done with 
model Q - might be too crude. 
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Fig. 5 Normalized distribution functions for j3, <fi and EW. 
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Fig. 6 Softness flux ratio against 2-10 keV flux for the 
AGN of our sample, where only an upper limit on the EW 
of the iron K a broad component is found (red squares) 
and AGN where a broad iron K a line is detected (yellow 
squares). 

In Fig.[6]we plot the 0.5-2 keV versus 2-10 keV soft- 
ness ratio for the objects of our sample, separating with dif- 
ferent colors objects where only an upper limit on the broad 
Fe K a line < 100 eV has been measured from objects where 
a broad line has been detected. The softness ratio should be 
sensitive to the presence of "warm absorbers", which im- 
print strong absorption features in the soft band 1 . We find no 
difference in the softness ratio in the two classes. We con- 
clude that it is unlikely that the number of iron K a broad- 
ened line detections is significantly affected by uncertainties 
in the underlying warm absorbed continuum determination 
(although in individual objects this effect may have an im- 
portant impact, as discussed by the above referred authors). 

4 Detection of relativistically broadened 
profiles in stacked spectra 

Relativistically broadened iron K a profiles have been de- 
tected in about 25% of the sources of our sample. In at least 
50% of the sample, the upper limits on the EW of a relativis- 
tically broadened profile are inconclusive as to establishing 
whether such a feature is present (cf. Fig.[3]and|4j. 

This statistical limitation, however, can be overcome, by 
stacking data together. In Fig.[7]we show stacked residuals 
for the objects of our sample, once they have been divided in 
four, equally populated (as of number of sources) 2-10 keV 
absorption-corrected luminosity classes. The stacked resid- 
uals were generated by summing together the residuals of 

1 This ratio is actually mostly sensitive to standard warm absorbers with 
photoionization parameters f; ~10— 100 and log(TVff) ~10 21-22 cm -2 
(Blustin et al. 2005). However, high-density, high-ionization warm ab- 
sorber are often found in conjunction with lower-density, low ionization 
counterparts. 
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Fig. 7 3-8 keV (3-6 keV in the insets) stacked residuals for the AGN of our sample, divided in four equally-populated 
2-10 keV intrinsic luminosity bins. The curves represent the: a) expected profile of an unresolved Gaussian profile, re- 
flecting the intrinsic pn camera energy resolution {dashed line); b) two maximally rotating Kerr profiles corresponding to 
EW=150 eV and 300 eV {dot-dashed profiles); c) two Schwarzschild profiles {long dashed-dot line) corresponding to the 
same values of EW. The luminosity bins are defined as follows: LI, Lx < 10 43 erg s , L2, 10 43 < Lx < 5x 10 43 erg s ; 
L3 , 5 x 10 43 < L x < 1.5 x 10 44 erg s" 1 ; L4, L x > 1.5 x 10 44 erg s" 1 



a fit with model Q, after the E^Gi and B((3, (f>, a) terms ter the profile is weaker. The integrated EW of the stacked 
had been removed. A clear broadened relativistic profile is profiles is <150 eV for any luminosity class (compare the 
present in the lowest luminosity sub-sample {Lx < 10 43 erg s~o)n Jrved profiles with the profiles of an extreme Kerr line in 
and possibly in the highest luminosity sub-sample. If the Fig.Qto Fig.[9}, therefore significantly lower than measured 
lowest luminosity bin is split into two equally-populated by Streblyanska et al. (2005) on the Lockman Hole AGN. 
bins, the broad line keeps being stronger in the lower-luminosity| 

bin (Fig ®. It seems to exist a sort of "Baldwin effect" 5 Can wfi acWeve ^ fm coy q{ a 

on the relativistic ally broadened component of the iron K a , . . ^ 

similar to what already claimed on the narrow component com P * * Samp e . 

of the same line (Iwasawa & Taniguchi 1993, Page et al. _ ,. , . . . _ , 

2003; see as well Jimenez-Bailon et al. 2005 for a different The P^ adl 8 m to e ^ la f n ° fthe 

. , , . in Active Galactic Nuclei (AGN) requires radiatively ef- 

mterpretation of the same data). _ . . • 

ficient accretion onto super-massive black holes (Lynden- 

Broadened profiles with comparable profiles are detected^ Bell 1969). In this scenario, the accretion flow reaches the 

in the stacked residuals of luminosity-matched samples of innermost regions close to the black hole, possibly the in- 

Seyfert Is and Seyfert 2s (see Fig. |9}, although in the lat- nermost circular stable orbit if it occurs in a standard geo- 
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Fig. 8 The same as Fig. [7] when the lowest-luminosity bin is split in two equally-populate (in source numbers) bins. The 
left panel corresponds to objects with Lx < 10 42 erg s -1 , the right panel to 10 42 < Lx < 10 43 erg s _1 . 
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Fig. 9 The same as Fig.0 for luminosity-matched sub-samples of Seyfert Is {left panel) and Seyfert 2s (right panel). 



metrically thin and optically thick disk. An inescapable pre- 
diction of this scenario is that relativistic broadening and 
skewing of emission lines should be common in AGN. The 
measure of the "true" fraction of objects where relativistic 
effects are detected would be crucial to establish if and how 
the standard AGN paradigm needs to be modified. Current 
AGN samples in X-rays are inevitably biased towards bright 
AGN, and furthermore towards bright Seyferts where the 
detection of a relativistically broadened iron K a profile was 
expected on the basis of observations prior to the Chandra 
and XMM-Newton launches. 

The AGN sample on which the study presented in this 
paper has been conducted is neither complete nor unbiased 
in any statistical sense. It is merely the collection of all the 
observations available in the XMM-Newton science archive 
at a given moment in time. This study needs to be extended 
(or restricted) to a statistically complete flux-limited sam- 



ple before any inference can be derived on the true fraction 
of AGN where relativistic effects are important, as well as 
on the properties of the innermost regions of the accretion 
flow around super-massive black holes (or, hopefully, on the 
black hole spin itself). The double bias outlined above needs 
to be removed. 

How can one build a complete sample, where the spec- 
trum of each member of the sample has enough statistics for 
meaningful constraints on the properties of relativistically 
broadened iron K Q lines to be derived? Let's consider again 
the EW versus net counts correlation in Fig.|2] For counts 
> 10 5 the correlation flattens. This represents the minimum 
number of net counts necessary to detect a line as strong 
as the average of the EW distribution. Moreover, a number 
of net counts >2x 10 5 ensures that a complete separation 
between measurements and upper limits is achieved. This 
condition is crucial to assess the significance of any non de- 
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tections. Hence, we conclude that in order to derive the true 
fraction of relativistic lines in a sample of AGN, each of the 
members of the sample has to be observed long enough to 
accumulate at least 200 000 net counts in the hard X-ray 
band. 

If we take the Piccinotti sample (Piccinotti et al. 1982) 
as example of a flux-limited sample in the hard X-ray band, 
only 6 out of its 27 type 1 objects have been optimally ex- 
posed with XMM-Newton. The detection fraction of rela- 
tivistically broadened iron K a lines in this sub-sample is 
50 ± 32%. We estimate that a full coverage of a hard X-rays 
flux limited sample such as the Piccinotti or the RXTE All 
Sky Survey (Revnivtsev et al. 2004) would require ~1 Ms at 
a flux limit ~2 x 10~ n erg cm~ 2 s _1 . This is large but not 
unreasonable time allocation for an XMM-Newton Large 
Program. 



6 Conclusions 

The main conclusions of this study can be summarized as 
follows: 

- relativistically broadened Fe K a lines are present in ~50% 
of well exposed XMM-Newton spectra 

- the detection rate and average properties are not signif- 
icantly different between type 1 and type 2 objects, al- 
though the analysis of stacked spectral residuals on lu- 
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minosity well-matched samples of Seyfert 1 s and Seyfert 2| 4 1 8 ; 927 



suggests a possible weaker line in the latter 

- the relativistically broadened Fe K a lines detected in our 
study cover a large range of EW. However, our stacked 
spectral residuals analysis finds typically EW<150 eV, 
far from the ~500 eV EWs claimed by Streblyanska et 
al. (2005) for the Lockman Hole AGN 

- relativistically broadened Fe K a lines are significantly 
more common in low luminosity (L x < 10 43 erg s _1 ) 
AGN 

- the average disk inclination is (cp) — 34°; the aver- 
age accretion flow power-law radial dependent index is 

</?) = -2.7 

- in 4 out of 6 objects where the measurement of the di- 
mensionless black hole spin yields meaningful constraints,! 
the Schwarzschild solution can be ruled out 

The conclusions of this study are still based on an in- 
complete and biased sample. However, a relatively mod- 
est investment of XMM-Newton time could allow them to 
be put on a firm statistical ground. At least 200 000 hard 
X-ray counts are necessary to ensure an unambiguous 
detection, as well as a clear discrimination between de- 
tections and upper limits. We estimate that around 1 Ms 
would be enough to optimally expose a hard X-ray flux limit 
sample at a flux limit ~2 x 10~ n erg cm~ 2 s _1 . 
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